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Abstract 

The two-photon two-neutrino interaction induced by magnetic field is investigated. In 
particular the processes 77 —* vv and 7 — > "fvv are studied in the presence of strong mag- 
netic field. An effective Lagrangian and partial amplitudes of the processes are presented. 
Neutrino emissivities due to the reactions 77 — > vv and 7 — ► jvv are calculated taking 
into account of the photon dispersion and large radiative corrections. A comparison of 
the results obtained with previous estimations and another inducing mechanisms of the 
processes under consideration is made. 
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1 Introduction 

Historically the reaction 77 — > vv was one of the first photon-neutrino processes considered 
in the context of its astrophysical application. In 1959 Pontecorvo suggested that (ev)(ev) 
coupling could induce reactions leading to energy loss in stars [1]. One of these processes, 



*E-mail addr ess : mch@ uni v . uniyar .ac.ru 
t E-mail address : mikheev@univ. uniyar .ac.ru 



77 - > vv, caused by this coupling was compared in [2] with other neutrino reactions and a 
rough estimation of the neutrino energy loss rate was obtained. In both papers the authors 
used the four-fermion (V-A) Fermi model. However, in 1961 it was proved that in this case 
the process under consideration is forbidden. This statement is also known as the Gell-Mann 
theorem [3] asserts that for massless neutrino and on-shell photons, in the local limit of weak 
interaction, the amplitude of the z/z/77-interaction is equal to zero. Really, in the center-of- 
mass frame the the left neutrino and right antineutrino carry out the total momentum unity in 
the local limit of the weak interaction. However, the system of two photons can't exist in the 
state with the total angular momentum equals to unity (the Landau- Yang theorem [4,5]). This 
statement also can be formulated in another way. The most general amplitude of the process 
can be written in the gauge invariant form: 

m = -% mi) T a „ v u t (-k 2 )] fffr, (i) 

vr V2 

where i is the neutrino flavour, i — e, //, r, and the tensors f a/3 = q a e lS — q^e a are the photon 
field tensors in the momentum space. The Gell-Mann theorem means that in the case when 
above- listed conditions are realised it is impossible to construct the nonzero tensor T a/3liu . Any 
deviation from the Gell-Mann theorem conditions leads to the nonzero amplitude of the process. 
For example, in the case of massive neutrino the process is allowed due to the change of the 
neutrino helicity [6,7]. The tensor T a p^ v has the following form in this case: 

Tafou = ( 2 $ie - !) -^Ih^^v (2) 

In the case of non-locality of the weak interaction the neutrino momenta, k 2fJi , become 
separated and the following structure arises [8-10]: 

Si { in? 3\ 1 

Tap,* = — ( In— y + - —5- [7afiV^0i - k 2 ) v + -f lx g V a{k 1 - k 2 )p] (1 +75). (3) 

It is seen that in both cases the amplitude is suppressed, either by small neutrino mass in the 
numerator or by large H^-boson mass in the denominator, and the contribution of this channel 
into the stellar energy-loss appear to be small. 

One more exotic case of nonzero amplitude is realised for off-shell photons [11,12], q^f^ 7^ 0, 
when the photon momenta can be included into the tensor T a p^ v 

i 1 

Tctpnv = ~— A ( 2 ^e - 1) — j 7 P (1 + 75) (epa,,,, qi/3 + q 2u ) (4) 
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However, there exist one more possibility to construct the nonzero amplitude (1). In the 
presence of the external electromagnetic field an additional tensor of electromagnetic field F^ v 
arises which allows one to construct the tensor T a j3^ v even in the case when Gell-Mann theorem 
conditions are realised. Therefore, one could say that external electromagnetic filed induces 
the two photons two neutrinos interaction. 

Regarding possible astrophysical application of the process discussed, it is interesting to 
consider the magnetic field case. It is well known that magnetic field presents practically in 
all astrophysical environments. In many of them the existence of intense magnetic field is 
assumed. For example, the typical magnetic filed of neutron stars is observed about 10 12 G. 
Some modern neutron star models consider the generation of magnetic fields up to 10 14 -r- 10 16 
G [13]. Moreover, the very recent observations of SGR and AXP pulsars indicate the existence 
of the magnetic field about 10 15 G [14]. Note that the strength of such magnetic fields exceeds 
essentially the so-called critical value B e = e/m 2 e ~ 4.41 • 10 13 G, which is a natural scale for the 
field strength 1 . It is known that strong magnetic field could enhance the processes suppressed 
in vacuum (see e.g. [15]), therefore it could be important to investigate the influence of strong 
external magnetic field on the process 77 — > vv. 

2 The process 77 — > vv in external magnetic field 

Previously the process under consideration was studied in the relevantly weak magnetic field, 
B <C B e . In the paper [16] an effective Lagrangian of the 777z/z/-interaction [17] was used to 
obtain the cross section and the emissivity of the process 77 — > vv with photon and neutrino 
energies much less than the electron mass. It was shown that the cross section of the process is 
enhanced by the factor (mv^/m e ) 4 (i?/i? e ) 2 in comparison with its counterpart in vacuum, where 
m\y and m e are W-boson and electron masses respectively. Another approach was developed 
in [18, 19], where an electron propagator expansion in powers of the magnetic field strength 
was applied to study the process 77 — > vv with energies greater than m e . In the low-energy 
limit the amplitude of the process obtained in [18,19] agrees with the result of Ref. [16]. In the 
paper [20] the results [16] and [18, 19] were slightly corrected. In particular, it was noted that 
the cross section of the process 77 — * vv has to be less by factor An. 

An investigation of the low energy two photon neutrino interaction in strong magnetic field 
was performed in [21]. The amplitude and emissivity of the reaction 77 — > vv was obtained 
in the four-fermion model without Z-boson contribution. 

The purpose of this paper is to investigate the two-photon two-neutrino processes in the 
presence of strong magnetic field with energies restricted only by the value of the magnetic field 

1 We use the natural units, c = h = k = 1, hereafter e > is an elementary charge. 
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strength, uj <C yeB. These processes are considered in the framework of the Standard Model 
using an effective local Lagrangian of the neutrino-electron interaction 

G 

where gy = ±1/2 + 2sin 2 9]y, 9a = ±1/2. Here the upper signs correspond to the electron 
neutrino {y = u e ) when both Z and W boson exchange takes part in a process. The low signs 
correspond to /i and r neutrinos (y = z/ M , z/ T ) when Z boson exchange is only presented in the 
Lagrangian (5), j a = v^( a {\ — 7 5 )z/ is the left neutrino current. 

In the third order of the perturbation theory the process 77 — > vv is described by two 
Feynman diagrams depicted in Fig. 1, where double lines imply that the influence of the external 
filed in the propagators of electrons is taken into account exactly. 



7(9") 




(7' 7") 



7(9') 



Figure 1: The Feynman diagram for the 1/1/77- interaction in magnetic field. 



The general form for the matrix element corresponding to diagrams in Fig. 1 is the following 

S = fI^J~ [ d 4 x d 4 y d 4 z Sp\(j'j)(qv — QAl5)S(x,y) x 

V2E'V 2E"V 2u'V 2u"V J lU,AJ 

x {e" 1 )S{y,z){e' 1 )S{z,x)}e- i ^' z -^ + (e',q> <-> e",q"), (6) 



where a ~ 1/137 is the fine-structure constant; q' a = (u/, q'), q" = (or,q' 
momenta of the initial photons with polarisation vectors e' a ,e'^ respectively; k a is the neutrino 
antineutrino pair four-momentum. S(x, y) is the electron propagator in the magnetic field 
which could be presented in the form 



" - are the four- 



S(x,y) = e^^S(x-y), 
y 

<f>{x,y) = -e d^ 



A^) + -F^-y)^ 



(7) 
(8) 
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where and are 4-potential and tensor of the uniform magnetic filed correspondingly. 
The translational invariant part S(x — y) has different representations. For our purpose it is 
convenient to take it in the following form 



S(X) = S-(X) + S+(X) + S ± (X), (9) 



where 



oo 

72 „ 





2 J2\ 



oo 

S ^ = -S/^7/(^ ( ' Y7W1 - ta,lh2T,x 



d 2 P = dpo^Ps, n ± = i(i±i 7l72 ), n| = n±, [n ± ,(A 7 )„] =o. 

Here 7a are the Dirac matrices in the standard representation, the four- vectors with the indices 
_L and || belong to the Euclidean (1, 2) subspace and the Minkowski (0, 3) subspace corre- 
spondingly, when the field B is directed along the third axis. Then for arbitrary 4-vectors A^, 
B p one has 

A*i = (0,^,^,0), < = (A ,0,0,A 3 ), 
(AB) ± = (AAB) = A 1 B 1 + A 2 B 2 , {AB\ = (AAB) = A B - A 3 B 3 , 

where the matrices A^ = ((pip)^, A^ = ((pip)^ are constructed with the dimensionless tensor 
of the external magnetic field, p>^ v = F^/B, and the dual one, (p^ v = \s^ upa (p p( j. Matrices A pv 
and A^ u are connected by the relation A pu — A pu = g pu = diag(l, —1,-1,-1), and play the role 
of the metric tensors in perpendicular (_L) and parallel (||) subspaces respectively. 

In spite of the translational and gauge noninvariance of the phase $(x, y) in the propagator 
(7), the total phase of three propagators in the loop of Fig.l is translational and gauge invariant 

y) + $(y, z) + x) = ~(z-x) ll F lu/ (x-y) v . 
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This fact allows one to define the amplitude of the process in the standard manner 

< ( 2xmwo (12) 

y/2E'V 2E"V 2uj'V 2uo"V 



where the amplitude M. can be presented in the following form 

Gf . a i ( -f-rV t-,A 



M = ^-j^ igy^y - g A U A \, (13) 



n y = e 3 



J d 4 Xd 4 YSp{ lfl S(X) lu S(-X - Y) lp S(Y)} x 
J d A Xd A YSp{ llxlb S{X) lv S(-X-Y) lp S(Y)} x 

x e ^XFY)/2 e iWX- q »Y) + ^ ^ ^ ^ (15) 



X 

„3 



with X = £ — rr, V = :r — y. 

In a general case substitution of the propagator (7) into the amplitude (13) leads to a 
very cumbersome expression in the form of the triple integral over the proper time. It is 
advantageous to use the asymptotic expression of the electron propagator for an analysis of the 
amplitude in strong magnetic field. This asymptotic could be derived from Eqs.(7)-(11) in the 
limit eB/\m 2 — p\ \ 1. In this case the parts S±, S ± of the propagator take the form 

§ (X) ~ — expf-^O / (p7)ll+m n e- 1 ^)" ( 16) 

b -^ A >- 2n 6XPl 4 >J (2tt) 2 V \ - m 2 ^ ' [ b) 

S + (X) ~ - JL [i( 7 9/9X)„ + m] <5 2 (X) n + exp(^) T(0, ^2), (17) 
^W--^^(X)^exp(-^i), (18) 

oo 

where T(a, z) is the incomplete gamma function F(a, z) — J t a_1 e _ *di. 

z 

By using, (16)-(18) the amplitude M. can be presented as a sum of the ten independent 
parts which can conditionally be divided into four groups: 1) S±S±S±; 2) S±S±S ± ; 3) S±S ± S ± ; 
4) S ± S ± S ± . Analysing these combinations one could expect that the leading on field strength 
part of the amplitude, namely ~ eB, arises from the combination S-S-S-. Really, substitution 
of S- in amplitude (13) gives (eB) 3 in numerator, while integration over transverse coordinates 
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X ± leads to the factor {eB) 2 in denominator. Then expending the amplitude in powers of 
inverse magnetic field strength B we obtain the linear on field dependence. However, it is easy 
to see that two parts of the amplitude (13) with photons exchange (e', q' <-> e", q") cancel each 
other exactly in this limit. Hence the amplitude of the process 77 —> vv doesn't depend on B 
in the strong magnetic field. 

The analysis shows that the independent on field contribution to the amplitude is given by 
the combinations S-S-S + , S-S-S ± and S-S ± S ± with all interchanges. One more contribu- 
tion comes from the expansion of the S-S-S- combination in the powers of inverse magnetic 
field strength B. Then substituting (16)-(18) into (13) we obtain the following result for the 
amplitude 

M ~ ^-j.e'y^gvU^ - g A n* up }, (19) 

n ^p = -^{(<?W) *w + v P » + 2 ((9" - 1') 1 )^ 

+ {q"l') p ^ - 1% (9V)„ + 1% {qip) P + /; p (q<p) v 

- l' vp {q"ip)»-I lxv {q"ip) p -I lxp {q'ip) v }, (20) 

ic^ 1 

+ W<pi') P w - {<j>i") P v (Mr + (<pi")tw (q<p) P + (sv)" 

- {(pl') vp (q"<p), - ((pi),, (q"<p) p - (<pl) w (<lW- (21) 
It is remarkable that the amplitude M. depends only on two types of integrals, I pu and n pup 

f d 2 p 

V = I»v(q) = -in J 7^pSp{7 /i 5||(p- ? )7^| l (p)}, (22) 
vr^p = -i?r y ■^^Sp{'j tl S ]] (p-q'')-f u S ]l { P )-f p S l] ip + q')}, (23) 

with 

s.M = #-±>_. 

p| — m 2 

Both types of integrals (22) and (23) can be presented in terms of analytical functions. The 
integral l pv can be written as: 



U?) = (v-^)#' 
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/4m|\ 
Id' 



where 

, . z 1 
H{z) = arctan —= — 1, z > 1, 

„, , If Z 1 + VI - -2 2 | 

w W = -2(7r^ ln wr^ + 2 -"7T^, l ' 2<1 ' 

The expression for 7r MJ , p can be presented in the following form: 

1 



nun 



n 2 n l2 n ll2 



(q'<pq"){(<pq)n(<pq")v(<pq') P n± 



Hw),(Aq"UAq') p H - {Aq)^q") v {Aq') p H" - (Aq) ,(Aq"U0q') p H'} 

Hq'q"UH^q"UH) P {H" - H>) + (qq")^q)^q"UAq') p (H - H") 

+ (qq')^q),(Aq'%^q') p (H' - H)] , (24) 

tt x = H' + H" + H 

n jX€ - 2rnl[qi(q>q"\H - qf(qq"\H> - q'^qq'\H"] 

qlq»qf-Aml[q»q'(>-{qiqi% ' 1 J 

where, e.g. iJ' = H{Am 2 e / q 1 ^). The result (19) can be also treated as an effective Lagrangian 
of photon-neutrino interaction in the momentum representation. Moreover it can be used to 
obtain the effective Lagrangians of axion-two photon (077) and three photon (777) interactions 
by making the substitutions 

\ A A A I n . ig a e x 

M ail =M{g v = 0, 9A^j=3i, -> 7^7- q») 

and 

G F 

.M 777 = M{g A = 0, gv^=-j^ -> e p ) 

V2e 

correspondingly. Here (? ae is the axion-electron coupling. 

Note also that the amplitude (19) in the low energy limit, uj <C m e , with gy — g A — 1 
coincides with the amplitude obtained in [21]. 

3 Contribution of the process 77 — > z/z/ into the neutrino 
emissivity of the photon gas 

To illustrate a possible application of the result obtained let us estimate the contribution of the 
process 77 — > vv into the neutrino emissivity of the photon gas in strong magnetic field. It is 
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convenient now to turn from the general amplitude (19) to the partial amplitudes corresponding 
to definite photon modes with polarisation vectors = (<pq) a / ' y/W, eL 1 ' = (<pq) a / s/o^ (i n 
Adler's notation [22]). These amplitudes can be written as 

2a Gp 1 

M« ± = -i 



tt ^2 

x{gv{[m'){m"), + (jq"U<iWW - 
tt V2 y/qfqf I 9jf v 

+ (W)ii[#vO)ii - g A (j^q)](H" - #')) - {wq')H"[g v {qq")\\ ~ 5a(?W)1 
-UWlH'[g v (qY) n +g A (q^q'')]\ (28) 



Then the neutrino emissivity (energy carried out by neutrinos from unit volume per unit 
time) can be defined as 

Q^ P = Qiiii + Q±ii + Qxx, (29) 
Qx'\" = (2tt) Va» ]T / |AW | 2 ^a» + £f ) 5 V + <z" - k' - k") 

d 3 q' , d 3 q" n d 3 k' d 3 k" 



X 



(2n) 3 2u> J v 7 (2tt) 3 2oj" J y ' (2tt) 3 2 ^ (2tt) 3 2 £f ' 



Here E[, E" are the energies of the neutrino and antineutrino of definite types % = v e , is^, v T \ 
uj', to" are the energies of the initial photons; f(ou) = [exp(uj/T) — l]^ 1 is the photon distribution 
function at the temperature T; the factor gyy = l — ^Syyi is inserted to account for the possible 
identity of the photons in the initial state. We would like to note that the integration over the 
phase space of initial photons in (30) has to be performed taking account for nontrivial photon 
dispersion law in the presence of strong magnetic field. Moreover, it is necessary to take into 
consideration the large radiative corrections in strong magnetic field which are reduced to the 
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wave-function renormalization factors Z\> and Z\" in Eq. (30). In the case of low temperature, 
T <C m e the dispersion law of the photon with polarization vector and corresponding 
wave-function renormalization factor Z ± could be presented in the following form 

where £ = ^ -jj- is factor characterising magnetic field influence. The photon with polarization 

vector Ea^ has almost vacuum dispersion law and wave-function renormalization factor in this 
limit, q 2 ~ 0, Z l} ~ 1. It is useful also to present here the element of the momentum space 

d 3 q = (1 + £)u 2 dujd(pdy, y = cos 9 y/l + f / ^1 + f cos 2 0, (32) 

where 0, are the polar and the azimutal angles. All of these facts lead to the rather complicated 
dependence of Q^ n ^ vV on magnetic field strength despite the amplitude (19) does not contain 
this dependence. In the case of the arbitrary parameter £ the value of the contribution into 
neutrino emissivity of the photon gas can be found only numerically, but when the strength of 
magnetic field is not too large, eB <C ml/a, one can obtain the analitical expressions for Qxy: 

Q = 8a2G } t 13 \g 2 v 7T 2 (7tt 2 C(5) + 330 C(7)) 

+ ^(9A(5) + 70vr 2 C(7) + 420C(9)) 

~ 8.5-10 7 T 9 13 ^^-, (33) 
y s • cm 3 



Qi 



127575 vr 3 mi 



g 2 v 24n 2 (l4vr 2 C(5) + 285 C(7)) 



+ g\ ((309tt 2 + 20)tt 2 C(5) + (1410tt 2 + 315)C(7) + 38430 C(9)) ] 



~ 6.9-10 8 T c 13 6rg 



9 q 

s • cm- 3 



(34) 



Qxx = 



8a 2 G 2 



il3 



127575 vr 3 m\ 



g 2 v n 2 ((609tt 2 - 50)C(5) + 10110 C(7)) 
10 



+ g\ ((879 7t 2 + 50)tt 2 C(5) + 5490tt 2 C(7) + 56700 C(9)) 



~ 3.3-10 8 T 9 13 



erg 
s • cm 3 



(35) 



Here T 9 is the temperature in units of V$K ~ 0.17m e , the effective constants g\ ~ 0.93, g\ ~ 
0.75 are summarized over all channels of the neutrino production, v e , v^, v T . The total neutrino 
emissivity in this case is 



Q 



LT 

r yy—*vu 



1.1 • 10 9 T, 



,13 er S 



s • cm 3 



(36) 



We also have made the numerical calculation of the neutrino emissivity caused by the process 
77 — > vv. In the low temperature limit T <C m e our result is represented in Fig. 2 where the 
emissivity Q^^ uP is depicted as a function of the parameter £. 




Figure 2: The low temperature, T m e , neutrino emissivity Q^^ uP dependence on 
the parameter £ = ^ -Jj- for different polarisation configurations of the initial photons; 
Qo = 10 8 T 9 13 erg/(s • cm 3 ). Short-dashed, dash-dotted and long-dashed curves correspond to 
Q\w\,Q±\\,Q±± respectively. Solid line depicts the total neutrino emissivity. 

The result presented by formular (36) and in Fig. 2 should be compared with the contribu- 
tions into the neutrino emissivity of the process 77 — > vv caused by the another mechanisms. 
For instance, the emissivity due to the finite neutrino mass is [7] 



Q^ vv ~ 1.4 • 10- 



: rpl\ 

1 9 



erg 
s • cm 3 



leV 



(37) 
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On the other hand, in the case of non- locality of the weak interaction, investigated in Ref. [10], 
one can estimate the emissivity, which is suppressed by the factor {m e /m w ) 4: : 

Q^^9.9-10- 10 T 9 13 ^. (38) 

S * CIX1 

It's obvious that the field-induced mechanism of the reaction 77 — > vv strongly dominates all 
the other mechanisms. 

It is interesting also to compare our result with the previous calculations of the neutrino 
emissivity due the process 77 — > vv in the weak and strong magnetic fields. Taking into 
account the remark by authors of [20] one could obtain from [16] the following estimation for 
the neutrino emissivity Q^^ uP in the weak magnetic field limit, B <C B e : 

Q^o,.io»r-(!) 2 ^. (39) 

It is evident that neutrino emissivity is substantially enhanced in strong magnetic field in 
comparison with its counterpart in the weak magnetic field case. 

In the limit B 3> B e the contribution of the process 77 — > vv into the neutrino emissivity 
was previously studied in [21], from which the following estimation could be obtained: 

Q*_ P ^0.7-10 8 T 9 13 ^. (40) 
" s • cur 

This result doesn't depend on the magnetic field strength B and it is at least ten times less 
than our result presented in Fig. 2. In our opinion, the authors [21] wrongfully didn't take into 
account photon dispersion and wave-function renormalisation in strong magnetic field. 

Let us note that in the presence of the magnetic field one more contribution into neutrino 
emissivity due to the process 7 — > ^yvv is possible. We would like to emphasize that it is 
the nontrivial dispersion law of a photon in the magnetic field that makes this reaction to be 
kinematically allowed. To our knowledge the process 7 — * 71/P has not been studied so far. 
Therefore it is interesting to compare the contributions into the neutrino emissivity from the 
77 — > vv and 7 — > '-/vv channels. 

To obtain the emissivity by the process 7 — > jvv one needs to make the replacements 
f{oj) — > (1 + /(^)) and q — > — q for one of the photons in (30). The analysis of the process 
kinematics shows that only one transition, 7y — > ^y±vv, gives the contribution into neutrino 
emissivity. The dependence of the neutrino emissivity due to the process 7 — > 71/P on the 
magnetic field strength is depicted in Fig.3. As is seen from both figures the contribution of 
the process 7 — > jvv into the neutrino emissivity turns out to be small in comparison with 
analogous contribution due to the reaction 77 — > vv in the case of not too strong magnetic 
field, B < 10 5 B e . 
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2 4 6 8 10 



Figure 3: The low temperature, T <ti m e , neutrino emissivity Q^^ lvV dependence on the 
parameter £. 

In the high temperature limit, T ^> m e , the main contribution into neutrino emissivity 
arises from the photons with parallel momenta squared are in the region m 2 e <C q 2 <C eB. In 
this case one can use vacuum dispersion law for both photon modes. We obtain the following 
estimation for neutrino emissivity in this limit: 



Let us compare this result with contribution into neutrino emissivity from process 7 — > vv. In 
strong magnetic field this reaction was studied in [23]. It was shown that the process 7 — > vv is 
kinematically allowed in the region q 2 > Am 2 . This lad to the suppression of the process in the 
low temperature limit by factor q 2 ~ Am 2 . In high temperature limit authors note that the main 
contribution into neutrino emissivity is defined by the vicinity of the point q 2 = Am 2 . However, 
the detailed analysis shows that the main contribution arises from the region m 2 <C q 2 (3. 
In this case the estimation of the neutrino emissivity coincides with the formula obtained by 
authors in the limit T 2 > eB,m 2 : 




(41) 




(42) 



Then the ration of the emissivities (41) and (42) is 




(43) 
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It is seen that the contribution into neutrino emissivity of the process 7 — > vv is always larger 
then 77 — > i/P reaction contribution, because -B ^> T 2 in strong magnetic field limit. For 
example, R ~ 5.86 for T = 2 MeV and B = 10(XB e . Nevertheless, we can see that these 
quantities have the comparable values. 

In summary, we have investigated the two-photon two-neutrino processes in the presence 
of strong magnetic field. The amplitude of the reaction 77 — > vv is obtained in a general case 
when the photons are not assumed to be on the mass shell. Therefore it can be treated as an 
effective Lagrangian of photon-neutrino interaction. We have also calculated the contribution 
into the neutrino emissivity due to the reactions 77 — > vv and 7 — * ^vv taking into account the 
photon dispersion and wave function renormalisation in strong magnetic field. We stress that 
in spite of the independence of the amplitude on the magnetic field strength B, the emissivity 
essentially depends on B. The comparison of our results with the other inducing mechanisms of 
the reaction 77 — > vv shows that strong magnetic field catalyses this process. Moreover it seems 
that the processes under consideration are the most dominating photon-neutrino reactions. 
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